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Abstract 

We present a derivation of the classical SIR model through a mean-field 
approximation from a discrete version of SIR. We then obtain a hyperbolic 
forward Kolmogorov equation, and show that its projected characteristics 
recover the standard SIR model. Moreover, we show that the long time limit 
of the evolution will be a Dirac measure. The exact position will depend on 
the well-know Rq parameter, and it will be supported on the corresponding 
stable SIR equilibrium. 



1. Introduction 

A very fruitful modeling paradigm in epidemiology is the so-called com- 
partmental models, with dynamics governed by mass-action laws. Most clas- 
sical epidemiological models are of this type, and this has led to a number 
of both quantitative and qualitative predictions in the disease dynamics 2]. 
More recently, there is a growing interest in discrete, agent-based, models 7]. 
See also [s| for a comparison between different models. In many cases, these 
models are thought to be more realistic, and able to capture important dy- 
namical features that are not present in the continuous models. 



Email addresses: chalub@fct.unl.pt (Fabio A. C. C. Chalub), msouza@mat.uff.br 

(Max O. Souza) 



December 16, 2009 



Here, we follow the ideas in [3], to study the the large-population regime 
of the discrete dynamics. In this way, we obtain a Hyperbolic Forward Kol- 
mogorov equation for the probability density evolution. It is well known, 
through the method of characteristics, that there is a strong linkage between 
solutions to first-order Partial Differential Equations (PDEs) and systems of 
Ordinary Differential Equations (ODEs). Thus, its not entirely surprising 
that projected characteristics of this PDE will be related to the classical SIR 
(Susceptible-Infected-Removed) ODE model. 

This modeling through PDEs has some advantages, in particular, it al- 
lows the introduction of higher-order effects, like, for example, stochasticity, 
adding a second-order term to the equation. On the other hand, if we con- 
sider a discrete model in population dynamics and consider its limit of large 
population (under suitable condition) we naturally obtain a PDE for p(t, x), 
the probability density to find the population at state x at time t. The ODE 
can then be obtained as the hyperbolic limit of the PDE, or, alternatively 
as the initial dynamics of the PDE. In short, this means that the dynam- 
ics of the discrete population can be approximated for short times and large 
population by a certain ODE — the derivation of this ODE requires an intro- 
duction of intermediate models, a stochastic differential equation or a partial 
differential one. 

The ODE approach can be seen in while the PDE modeling was the 
subject of a previous work from the authors, where the replicator equation 
was obtained as the limit of the finite-population discrete Moran process [i[ . 
The resultingequation is of singular type and required a specific analysis of 
its behavior [4|]. 

In this work, we will study in a certain level of detail the SIR epidemic 
model from the PDE point of vie This is one of the most elementary and 
well studied model in mathematical epidemiology. See 0,0]. 

In particular, we shall prove that the solution of the SIR-hyperbolic- 
PDE obtained as a first order expansion in the inverse of the population 
size from the discrete-SIR converges when t — > oo to a Dirac-delta measure 
supported at the unique stable equilibrium of the SIR-ODE. Moreover (and 
this will be proved in a forthcoming work) the solution of the SIR-hyperbolic 



lr The expression SIR appears in the literature in two different context: one as a discrete 
evolutionary system, used in general in computer simulations; the second as an ODE 
system. We hope that all these different meanings are clear from the context. 
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is approximated by the SIR-ODE in all time scales and approximate the 
discrete-SIR for short time scales. This provides a framework to unify all 
these descriptions. In a forthcoming work, we will also introduce the SIR- 
parabolic-PDE which has the inverse behavior (approximate the discrete-SIR 
for all time scales and the ODE-SIR for short times). 



2. Discrete and continuous SIR models 

Consider a discrete SIR model, i.e., consider a fixed size population of 
N individuals, each one in one of the three states: n individuals Susceptible 
(i.e, individuals with no imunity), m individuals Infected (i.e., individual 
currenctly infected by a given infectious disease and able to transmit it to 
the susceptibles) and N — n — m Removed (after infection idividuals have a 
temporary imunity and then are removed from the dynamics; after certain 
time they become susceptible again). This is a very simple model for non 
lethal diseases transmitted by contact, e.g., normal influenza. At each time 
step of size At > we select one individual at random: 

• If it is S, then it changes to I with probability proportional to the 
fraction of I in the remainder, am/(N — 1); 

• If it is I, then it changes to R with constant probability j3; 

• If it is R, then it changes to S with constant probability 7. 
This can be summarized in the following diagram: 

S + I^I + I, I-^R, R^S. 

This model is also called SIRS, and when 7 = we recover the classical 
SIR model. For simplicity, however, we will call it the SIR model in this 
work and all results presented here include the case 7 = 0. 

Constants a, (3 and 7 depend, in principle, in N and At. As we are 
interested in the limit behavior when N — ► 00, At — > we will assume the 
following scaling relations 

lim — — = a , lim — = b , lim 



AT-»oo,At->0 NAt iV->oo,At->0 NAt iV-»oo,At->0 NAt 

with ab ^ 0. For further informations on scalings, see 
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Let -P(jv,At)(£j n i m ) be the probability that at time t we have n susceptible, 
m infected and N — n — m removed, where the total population N is constant 
and the time step is given by At > 0. Therefore 

( n _|_ l)(m — 1) 

P(N,At) (* + At, n, m) = a — — - - — -P(jv,At) (t, n + 1, m - 1) 

777/ H - 1 TV — Th — 771 -\- 1 



+ 



N 



iV- 1 



m . ^ x N — n — m , 
+ Tr(l-P) + ^ 1-7 



N 



N 



Now, define x = n/N, y = m/N and p(t, x, y) = P(t, xN, yN; N). Then, 
using p(t, x,y) = p and keeping terms until order 1/N: 



pit + At, x, y) = a 



fc+j) (y--k) 



AT 



A/ 



A^y " V 7 AT' 



+ (* (i - j^r) + 2/(1 - 0) + (i 



x-y)(l-7) }p(t,x,y) 



p + — [(a (y - x) + f3 + 7) p + (axy - 7(1 - x - y)) d x p + (f3y - axy) d y p] 
P+jf[d x ((axy - 7(1 - x - + d y (((3 - ax)yp)} 



Finally, 

d t p = d x ((axy - c(l - x - y))p) + d y ((b - ax)yp) . 
subject to probability conservation, i.e, 



(1) 



dt 



p(t, x)dx = 0. 



(2) 



3. The SIR model as a transport problem 

Let x = (x,y), and let $t(x) be the flow map associated to the SIR 
system 

X = c(l - X - Y) - aXY , 
Y = (aX- b)Y , 
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then, if p (x) is C^R 2 ) function, and p(t,x,y) = e Q W- Q ^-^p {<$>_ t {x)), 
with x = (x,y), and Q satisfies 



where F denotes the right hand side of the SIR system. Let S denote the 
unit simplex in R 2 . 

Fix x £ S. Then, we have that 



= -F(0 t (x o )) • VQ(<t> t (xo)p(t, <f> t (xo)) + d tP (t, t (x o ) + F(0 t (x o ))Vp(t, t (x o )) 
= V ■ Fp + FVp + d t p 
= d t p + V ■ (pF) . 

Thus, the SIR system are the characteristics of ([T]), and the probability 
density should be transported along them. We now make this calculation 
more precise. 

Definition 1. Let F : U C lR n — > M. n , be a Lipschitz vector field, where 
U is an open set, and let Q C U be compact. We say that Q is regularly 
attracting for F, if there is an open set V with a piecewise smooth boundary 
and Q C V C V C U, such that, if $ t denotes the flow by F restricted to V , 
then we have that uj(V) C 0. 

Theorem 1. Let Q be a domain with piecewise smooth boundary dQ. Let 
F : U C M. N — > M. N be Lipschitz, with Q C U being a regularly attracting set 
for F. Let pq G L x (Q) be nonnegative. Then the equation 



F-VQ 



V-F, 



e -Q(^(xo)) p(t;0t(xo))=e Q(xo) po(xo) , 



Hence, 



= e Q(^(x )) d [ e -«(^( x o)) p ( t; t ( Xo ))] 



d t p + V • (pF) = 0, p(0, x) = po(x) 
has a unique solution 



(3) 



p(t,x) = e Q(x) - Q( *- t(x)) po($-t(x)) 

Moreover, p is nonnegative, and supp(p(t, •)) C Q. 
In addition, z/supp(po) C fi, then 




(4) 
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Proof. Existence can be shown as follows by considering the weak formula- 
tion, Let W = [0, oo) x ft, and let ip G C C (W). 

// p(t,x)9t^(t,x)dxdt + // p(t,x)V^(t,x)dxdt + 
J Jw J Jw 

+ / p(0,x)^(0,x)dxdt = 0. (5) 
Jq 

Choose ip G C c ((0, oo) x ft). Then fl^D becomes 



p(t,x)d t ip(t,x)dxdt + // p(t,x)V^(t,x)dxdt =0. 

Let y = <p t (x) and 7ft (cc) = det(d x $ t (x))- Also let PVf = $ t (W). Then, the 
first integral becomes 

f f e Q ^y»- Q ^p (y)dMt, $ t (y))r) t (y)dydt. 
The second integral becomes 

/ / e c 3^^)"«^ Po (y) v ■ ^(*, $ t (y))^(y)dydt. 

Combining both integrals, we can write 

/X eQ( * t(y)) ~ Q(y Wy)^(*, $*(y))^(y)dydt. 

On integrating by parts, we have that 

I f e -Q(y) po (y)ip(t, $ t (?/))e Q( * t(y)) [F ■ VQ - V • F] dydt = 0. 
J Jw t 

We have that p is clearly nonnegative. Moreover, Let V = supp(u(t, ■)) and 
let Vt = (&-t(W)). Let ip be a vanishing function in ft. Then 

0= / M(t,x)^(x)dx= f e Q(<I>t(y) - Q(y) Mo(y)^(y)div(F)dy 

Hence supp(w(t, •)) C $t(ft). Since ft is regularly attracting for F, we have 
that $t(ft) C ft. If supp(po(y)) C ft, then we extend po by defining it to be 
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zero in M. — U. Then we have that 

s/ n M(t ' )dx = sX M(v)dx 

= - / div(u(i, -)F)dx = 0. 
Jv 

□ 

Let 

F s (x, y) = (c(l - x - y) - axy, axy - 6y). 
We define the reflected SIR field by 

F RS (x,y) = (c(l - x - y) - axy, (ax - b)\y\) 

We immediately have 

Lemma 1. Let Q be the simplex in the nonnegative orthant of R 2 . Then 
Frs is a C 1 vector field in M 2 , and Q is regularly attracting for Fr$- 

Thus for smooth solutions, we have 

Theorem 2. Consider the Cauchy problem for (CQ), with a L 1 non-negative 
initial condition po . Then there exists a unique solution satisfying (j2J) • More- 
over, p{t, x) > 0. 

4. Asymptotic Behavior and measure solutions 

We recall that the dynamics of SIR are controlled by the parameter R : = 
a/b (see i.e, 

Proposition 1. Let xi = (1,0), and X2 = yjj^, ^ ^1 — ito)) ^ e ^ e ^ wo 
equilibria of SIR. xi referred to as the disease free equilibrium and X2 as 
the endemic equilibrium. If Rq < 1, then any solution that starts in the 
nonnegative orthant ofW 2 approaches x for large time. If Rq > 1, then Xi 
is the limiting point. 

With this point of view we have 
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Theorem 3. Let p be a solution of (JT]) ; satisfying (j2J). Then, in the Wasser- 
stein metric, we have that 

5v, , Rn < 1 



t— >oo^^ ' ^ 1 5 X2 , -Rf) ^> 1. 

Proof We deal with the case Rq < 1; the case i? > 1 is analogous. Since 
-Ro < 1> we have that xi is the globally asymptotic stable equilibrium. Then 
given, 5 > 0, we can find T > 0, such that, for t > T, we have that 

$t(S) C B,( Xl ). 

Let ^>(x) be a continuous function. Then, for t > T, we have that 
/ p(t, x)^(x)dx = / x)^(x)dx. 

JS Jb s (xi) 

But, let e > be given. Since tp is continuous, we have 5 > 0, such that 

^(Xi) - e < / p(t, x)^(x)dx < ^(x ) + e, 
and this proves the claim, since we have that 

/ p(t, $ t (y)dy= / e Q (* t ^)" Q W^ (y)^(^(y))^(y)dy 



□ 
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